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Abstract 

This study deals with the further development of nuclear spin model of scalable quantum register, which 
presents the one-dimensional chain of the magnetic atoms with nuclear spins 1/2, substituting the basic 
atoms in the plate of nuclear spin- free easy-axis 3D antiferromagnet. The decoherence rates of one qubit 
state and entanglement state of two removed qubits and longitudinal relaxation rates are caused by the 
interaction of nuclear spins-qubits with virtual spin waves in antiferromagnet ground state were calculated. 
It was considered also one qubit adiabatic decoherence, is caused by the interaction of nuclear spin of 
quantum register with nuclear spins of randomly distributed isotopes, substituting the basic nuclear spin- 
free isotopes of antiferromagnet. We have considered finally encoded DFS (Decoherence- Free Subspaces) 
logical qubits are constructed on clusters of the four-physical qubits, given by the two states with zero total 
angular momentum. 

Keywords: Easy-axis antiferromagnet, decoherence, indirect coupling, inhomogeneous magnetic field, nu- 
clear spin, quantum register and qubit. 
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1 Introduction 

In early papers [1-4], we have considered a model of NMR quantum register, which is based on the nuclear 
spin-free easy-axis 3D antiferromagnet at low temperature in homogeneous field. It was shown that the range 
of indirect coupling can ran up to a great value close to critical point of spin-flop quantum phase transition in 
antiferromagnet. We have extended the previous model into the case of inhomogeneous external magnetic field. 

It was proposed to use the natural antiferromagnetic crystals with easy-axis anisotropy as an antiferromagnet 
thin plate (or film). As examples, they may be crystals CeC2 with tetragonal and FeCOa (siderite) with trigonal 
symmetry. The basic isotopes of these crystals ^^C, '^^fe (91.7%), ^^O, ^^°'^''^Ce (99.6%) have no nuclear spins 
(in brackets percent isotopic abundance is given). To form the one-dimension nuclear spin chains the isotopic 
substitution atoms, such as ^^C in corresponding crystal lattice sites, for isotopes ^^C with nuclear spins 1/2, 
are proposed. One would expect that period of such solid state NMR quantum registers may be much more 
than periods of crystal lattice. 

The simple antiferromagnet model to be studied here consists of two incorporated to each other tetragonal 
magnetic sublattices A and B with N — N±Nz sites in each sublattice, where N± = N^Ny 3> 1 are the sites 
numbers in plane of plate (x,?/-axes) and > 1 is the sites numbers in z direction. The atom sites of sublattices 
are numbered respectively by numbers j and i. Each sublattice constant in the plane of the plate is a± and 
along symmetry axis is a^. 

The external magnetic field B (x) in considered model (Fig. 1) is directed parallel to z-axis and to the 
antiferromagnet easy axis. It is assumed here, that the field gradient is of the order dBz (x) /dx = G ~ 0.1 T//im 
(x— axis in plane of the plate along nuclear spin chain). This value of the gradient corresponds to the difference 
of resonance frequencies of the order of 100 kHz for two nuclear spins, being separated by 100 a± (aj_ ~ 1 nm). 

The qubits number in quantum register will be limited by planar structure dimensions. For example, the 
structure with linear dimension of the order of 10 fim will have 100 qubit register with period L = lOOaj^. 
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Figure 1: The scheme of antiferromagnet based nuclear spin quantum register in external field lower than critical 
field for spin- flop phase transition B (x) <C Be- The oriented (that is nonprecessing) arrows represent here the 
ground states of corresponding individual Bloch vectors. The nuclear spins 1a (shot red arrows) are contained 
here only in atoms of sublattice A. 

The starting spin Hamiltonian of 3D easy-axis antiferromagnet with interaction only between neighbouring 
atoms, which belong to the distinct sublattice, is represented for our model as 

(N ^ \ 

J2 B (x,) Saz {Ti) +Y,B (xj) Sbz {rj) j + (1) 

N Z 

+2-fsh/Z ^ ^ { 5a (r^) Sb in + rs) + Ba Saz (r.) Sbz (r^ + rs)}, 

i 5 

where Sa (ri) and Sb {ti + rg) are electron spin operators {S = 1/2) for neighbouring sites of sublattices A 
and B, B {xi) — B + Gxi, B is the field value at the origin of the coordinates Xi, Z = 6 is the number of 
neighbouring atoms for tetragonal sublattice, 75 = 175.88 rad GHz/T (75 /27r = 28 GHz/T) is electron spin 
gyromagnetic ratio, Ti — 1.054 • 10"'^^ Js/rad. 

The direct product of spin operators in Eq.(l) written in matrix representation is common designated by 
symbol (8). In the following this symbol will be for brevity omitted. 

In Eq.(l) the parameters Be ^ (10 - 100) T, Ba ^ (lO^^ - 1) T > 0, Be = ^2BeBa + B\ are ex- 
change field, anisotropy and critical spin- flop field for easy-axis antiferromagnet (particularly, for FeCOs: 
Si3 = 35T, = 3.3T ft! 15.5T). 

It was shown in Ref. [1-4] , that indirect interaction of two nuclear spins in the model of antiferromagnet-based 
nuclear spin quantum register with inhomogeneous magnetic field essentially grows and qualitatively changes its 
character, if the value of local field in the mid point of two considered nuclear spin is close to the critical field for 
quantum phase transition of spin- flop type in bulk easy-axis antiferromagnet. The corresponding coordinates 
of nuclear spins were named by us as turning points. 

In present paper, we have presented the further investigations and development of this model. With the more 
refined results of the model analysis it was investigated here in details the one-qubit and two-qubit nonadiabatic 
decoherence and longitudinal relaxation rates are caused by the interaction of nuclear spins with virtual spin 
waves in antiferromagnet ground state. It was considered also one qubit adiabatic decoherence, is caused by the 
interaction with nuclear spins of randomly distributed isotopes, substituting the nuclear spin-free isotopes of 
basic antiferromagnet. We have considered finally, as an example, encoded DFS (Decoherence-Free Subspaces)- 
states of logical qubits are constructed on clusters of the four-physical qubits, given by the two states with zero 
total angular momentum. 

The paper is organized as follows. In Section 2 we have considered some general expressions of one qubit 
decoherence and relaxation processes are caused by interaction of nuclear spin with spin waves in easy axis 
antiferromagnet plate with inhomogeneous magnetic field. In Section 3 it was considered the nonadiabatic one 
qubit decoherence and longitudinal relaxation. In Section 4 it was considered the decoherence of two qubit 
entangled quantum states. In Section 5 we have considered the adiabatic decoherence is caused by interaction 
with nuclear spins of random distributed isotopes, substituting the spin-free atoms of basic antiferromagnet. In 
Section 6 we have discussed an encoded DFS (Decoherence-Free Subspaces)-states of logical qubits on clusters 
of the four-physical qubits. In Conclusion we discuss the some prospects of considered quantum register model. 
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2 Some general expressions for one qubit decoherence and longitu- 
dinal relaxation in antiferromagnet-based quantum register 

The antiferromagnet electron spin system of considered model plays the role of an environment for nuclear spin 

quantum register, whose interaction with spin wave leads on the one hand, to indirect coupling between nuclear 
spins and on the other hand to decoherence and relaxation processes of their states. 

Let us consider the decoherence and relaxation processes of quantum state for single nuclear spin placed at 
position k on axes x of sublattice A which is caused by its interaction with virtual magnon excitation in anti- 
ferromagnet. Such processes are described by transverse and longitudinal relative to external field components 
of Bloch vector (Ref.[5]) 

P- (fc, r) = (P+ (fc, r)) * = (fc, r) - iPy (fc, r) = 2 Tr (/" (fc) p (fc, r)) = 2 Tr/ (/" (fc) pi (k, r)) 

{k, r) = 2 Tr, (/, (fc) Pi {k, r)) , Trj if = 1/2, Tr, 7+7" = 1, (2) 
where non-steady reduced to one nuclear spin density matrix in antiferromagnet is represented by 

PI (fc, r) = 1/2 [1 + P, (fc, r) h {k, r) + P- (fc, r) /+ {k, r) + P+ (k, r) 7" (k, r)] (3) 

Let us pass next in Eq. (1) to dimcnsionless designations: 

hs = Hs/hujE, IsBe =coe,0< Ba/Be = 6^ < 1, < B/Be = b, Ga±/BE =g^ 10"^ . (4) 

The interaction of fc-th nuclear spin with external magnetic field and magnon excitations will be described 
here by dimcnsionless Hamiltonian with hyperfine interaction of the form (Ref.[4]) 

h = H/hcoE = hs- {oJi (fc) - aS^ {k)) h {k) + a/2 (7+ (k) S' {k) + 7" (fc) S+ (fc)) = 

= hs~ {i^i (k) - a (1/2 - ^) 7, (fc) + Ahis (k)) , (5) 

where 7^, 7=*= = 7^, ± ily, Sz, = Sx i iSy are nuclear and electron spin operators, ui (k) = {ji/^s) bk is 
resonance nuclear frequency in local field bk = b + gk, 77/75 ^ 10~^, a = A/coe ^ 10~^ — 10~* is isotropic 
dimcnsionless constant of hyperfine interaction and 

(0|l/2 - 5, (fc) |0) « (0|5+ (k) S- (k) |0) = ^ (fc) < 1. (6) 

is known as "spin contraction" . 

The perturbation Hamiltonian, corresponding to the relaxation and decoherence processes in isotopic pure 
antiferromagnet, has the form 

Ahis (fc) = A/^S {k) + Ahf^ {k) , (7) 

where 

a4s^ (fc) = a/2 (7+ (fc) S- (fc) + 7- {k) S+ (k)) (8) 

The second term in Eq.(7) 

Ahf^ (fc) = ah {k) {S, (fc) - (0|5, (fc) |0)) « -al, {k) ( {k) S+ (k) - ^) (9) 

describes two-magnon interaction, which is similar to two-phonon interaction (Ref.[5],§3.4). This mechanism 
causes in particular the modulation of nuclear spin resonance frequency without changing its state (adia- 
batic decoherence). It leads to the temperature depending decoherence rate, which is negligible small at 
(^bc — b^ hysBc ) ' ^^^^ ^he system is close to ground state. Therefore, we will neglect next the contribution of 

terms Ahfg (k) and will use as the perturbation Hamiltonian only the expression Ahis (k) = Ahfg {k). In this 
case, relaxation of transverse component of Bloch vector is accompanied by nuclear spin flopping (nonadiabatic 

decoherence). At the same time, the relaxation of longitudinal component of Bloch vector also occurs. Thus, 
these two processes may be considered here as one unified process of nuclear quantum state damping. 
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We assume now that interaction of nuclear spin, which is initiahy at coherent state (with nonzero nondiagonal 
elements of density matrix p/ (fc, 0)), and antiferromagnet in ground state, is turning on at the initial moment r = 
0, when nonperturbed density matrix is represented as direct product p [k, 0) = pi [k, 0) ps (0) = pi {k, 0) |0)(0|. 

Let us go next to interaction representation for density matrix relatively to Hamiltonian ho (fc) = hs — 
{ui {k) - a/2) , 4 {k): 

Pin [k, t) = exp {iho {k) r) p {k, r) exp {-iho (fc) r) (10) 
and to the equation for density matrix of nucleus-electron system 

idpin {k, t) /dr = [A/i/s {k, t) , p„ (fc, r)] , (11) 

where 

A his {k, t) = exp [iho [k) r] A/i/s [k) exp [-iho (k) r] = (12) 

= a/2 exp {-i {wi (fc) - a/2) r) 7+ (fc) S' {k, r) + H.c. 
It is follows from Eq.(ll) in the second order of perturbation theory 

T 

idpin {k, t) /dr « [Ahis {k, r) , p {k, 0)] - i J [Aft/s (r, k) , [Ahis {k, r') , p (fc, 0)]]dr'. (13) 



Let us write the derivative of expression (2) with respect to time, by using Eq.(13) and perform then the cyclic 
permutation under tracing. Finally, accounting the relation 

Tr ([/- (fc) , Ahis {k, r)] ps (0)) = 0, we will find 

d [P- {k, t) exp (i {uJi {k) - a (1/2)) r)] /dr = 2 Tr{r (fc) dpin {k, r) /Or} « (14) 

T 

^-2Tv{J [[l-{k),Ahis{k,T)], Ahis{k,T')]p{k,0)}dT'. 



By defining the transverse Bloch vector component in the form 

P- {k, t) = P- (fc, 0) exp [-i [wi {k) - a/2) r - 7_l (fc, r)] , (15) 

we will have: 

d [P- (fc, r) exp {i {uJi (k) - a/2) r)] /dr = (16) 

= —d'y± {k,T) /dr ■ P~ (fc,T)exp(i {coi (fc) — a/2)T) , 

where Re 7^ (r) is dccohcrcncc decrement and Im7j^ (t) is phase shift. 

We will represent the nuclear density matrix in the right part of Eq.(14) by expression 

PI {k, 0) « 1/2{1 + 2P, {k, 0) h (k) + P- {k, 0) 1+ (k) + P+ {k, 0) /- (fc)}|0)(0|. (17) 

Using in Eq.(14) the perturbation Hamiltonian (8), ignoring the factors 
exp{±i{uji (k) —a/2)T) and taking in the context of second order of perturbation theory in the left-hand 
side of Eq.(16) P^ {k, r) « P~ {k, 0), for the decoherence rate we will obtain 

r 

dRej^{k,T)/dT = ReTri J (0| [ [/" (fc) , A/ijs (fc, r)] , Ahis {k,T')] 1+ {k)\0)dT' = (18) 



= y 2 Re j {0\S- {k, t) S' {k, r') + S+ {k, t') S' {k, r) |0)rfr'. 



Let us consider next the relaxation longitudinal component of Bloch vector. We will write: 

dP, {k, t) /dT = 2Tr h^dpin {k, t) /dr « 



(19) 



T 

«-Re2Tr J [[h,, Ahis ik,T)] , A/i/s (A;, r')] (A, 0)|0)(0|dT'. 



Taking {k, r) = Pz {k, 0) exp (— 7|| {k, r)) , in the second order of the perturbation theory we will obtain: 

dPz {k, r) /dr « -d7|| (fc, r) /dr ■ P, (fc, 0) , (20) 
where the longitudinal relaxation rate is 

(i7|| (A:, r)/dr= (dill (A:, r)/dr)*= (21) 



T 

2ReTiiJ (0| [[h {k) , Ahis {k, r)] , Ahis {k, r')] h (k) \0)dT'. 



As a result of Eq.(21) we will have 



{k,T) jdi 



(22) 





whence it follows that for the considered mechanisms the rate of relaxation of longitudinal component is equal 
to the rate of relaxation of transverse component (decoherence rate). 

3 Nonadiabatic decoherence and longitudinal relaxation rates of one 
qubit quantum states 

For the calculation of nonadiabatic one qubit decoherence and longitudinal relaxation rates we used next the 
results of antiferromagnet Hamiltonian diagonalization, obtained in Ref.[4]. The transverse components of 
electron spin operators will take the form 



S- {k,T)={S+{k,T))^ = 



(23) 



= J [u* {q±, E) exp {iE_T) ^+ [qy, E_) + v {q±, E) exp {-iE+t) ^ {qy, £+)] exp {iq^k) dEdq±, 

where {qy,EJ), £^{qy,E^) are operators of creation and annihilation of spin magnons for two branches of 
magnon states, which propagate along direction .x-axis with dimensionless energies E± = E±h, Eisa, continuous 
energy parameter, and wave vector component qy in the range of qy to qy+dqy. The transformation coefficients 
have here the following asymptotic form (Ref.[4]) 



u* {q±,E) = 



V* {q±,E) 



E 



+ 1 • exp 



E 



exp 



Eqx - j E{q±) dqx 



Eqx - / E{q±) dq^ 



0{9).. 



0{g), 



(24) 
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where be = Bc/Ba is dimensionless critical field of phase transition and E {q±_) + 

Let us restrict next to low magnon excitation mode with energy Further, let us insert week magnon 
damping s {E- — )■ E- + is, E- ^ s > 0) and make a set of rearrangements, after which we will obtain 

Red7_L(fc,T)/dTW (25) 



Re 2 / j u* (g_L, E) u {q'^, Qy, E) exp (- {iE- — s) (r' - r)) exp (z {q^ — q'^) k) dq'^dq±dEdT' . 



For the estimation of magnon damping it may be used the line width of antiferromagnetic resonance s ~ 
js^B/ujE. The typical value of APR line width is AB ~ 10"^ T and s ~ AB/B ~ W'^. 
Upon integrating over r', we have 

Red'jj_{k,T) /dr K (26) 



/rT2^®2 / (9-L,-E')w(9x>9j/>-E')exp(«(gx - g^)''^) 
(47r) J 



, 1 - exp {{iE- - s) r) 



-iE_ + s 



dEdq±dq'^. 



By using the expression E (q±) « y 6^ + we write 

dq± = 2nq±dq± = 247r£' {q±) dE {q±) , < g_L < tt 



(27) 



and take notations Ab^ = be — b — gk, ^ = E [q^^ — be- Upon integrating over E and q'^ (Ref.[4]), we will 
transform Eq.(26) for decoherence rate to the form 



3a^ 

Red7_L (fc,T) /dT = —R± (A6fe,r) 
Ztt 



(28) 



2 V^+^-bc 

^ y (^1 + bl + bc + ^^ Y (^, +A6fe, r) dC, 



where 



Y{^ + Abk,T) = Re 



1 - exp [{i ($ + A6fc) - s) t] 
+ Abk) - s) 



(29) 



(^ + Abk) sin ((C + A6fc) r) exp (-sr) 1 - cos ((^ + Abk) t) exp (-sr) 



(e + Afefe)^ + s2 



(e + A6fc)^ + s2 



For the case of Abk s > 0, after omitting in dominator the Eq.(28) take the following explicit form. 



{Abk,T)^{{ Jl + bl + bc- Ab 



si (A&feT) - si I I + — - 6c + A6fc I r 



+ 1/7 



COS (A6fer) - cos ( | \/ ^'c + ~ + Afe^ j r 



} exp (— sr) 



1 + 62, + 6c - A6fc { 



be 



Abk 



bl + ^-bc + Abk 



cos (A6fer) 
A6l 



cos 



6c + A6fc r 



+ r I si(A6feT) - si ( ( \/6^ + — - 6c + A6fc I r 



^62, + ^ - 6c + Abk 
exp(-ST)}+ (30) 



+ 



ci (A6feT) - ci I ( \/^c + ^ -hc + Abk 



exp (— sr)} > 0, 



where si (a;) = — J ^^dt, ci{x) = — J ^^dt are sine-integral, and cosine-integral, and 

X X 

lim [ci (ar) — ci (6t)] = log j. 

T— 7-0 

Fig. 2 gives the r-dcpcndcnce of R± (Abk,T) for distinct values of parameter Abk- 

We notice, that oscillating part of the rate of decoherence tends fast (microseconds) to constant. 

Taking then r^oo(r>'l/s'~ 10^), we will obtain for inverse time of decoherence 

3^2 

1/Td = WBRed7_L (A:,oo) /dr = wb-— i?_L (6fe,oo) = (31) 

ZTT 

= '^SSl^ { , + log TT }• 

Abk^bl + ^-bc + Abk 

For values ~ 10~^ - 10~^, s ~ 10~^ < A6fc «; Ifc, WE/27r ~ 10"Hz decoherence time is Tjj ~ 
Abk/iOEsa^ sec. 

It follows that decoherence time, caused by one-magnon processes near turning points tends fast to the 
low value (for Abk > 10~^ decoherence time may be from milliseconds to seconds). Note that character of 
nonadiabatic decoherence rate depends on the anisotropy of antiferromagnet (through parameter be) and from 
inhomogeneity of external field (trough parameter Abk)- 

The expression for frequency shift is 

Im rf7_L {k, oo) /dr = 



2 



^Jbl + ^-bc + Abk 



= -^{(Vl + ^^ + ^c-AUlog V - - - ■- ^yjhl + --bc}^ (32) 



«-^(Vl + ^c + &cjlog^^^ 
It does not depend from magnon damping parameter. 



4 The decoherence of two qubit entangled quantum states 

The arbitrary state of pair spin-qubits in quantum register with zero Bloch vector values P (k, t) =2 Tr/ I^pi {I, k, t) 
P {I, r) = is described by the following reduced density matrix of nuclear spin system {a, = x,y, z) (Ref.[5], 
§§2.5-2.7): 

PI (/, k, t) = Tri,...,m^i,k,...N PI (1, I, ...k, ...AT, r) = 1/4{1 + ^ 40^,0 {I, k, r) (/;,„ 4^)}. (33) 

The non-steady pair entanglement state will be determined in the form 

PI {I, k, t) = 1/4 {1 + {I, k, t) {IiM + G+'- {I, k, t) (Ifl^) + G-'+ {I, k, t) (/;+/,-)}, (34) 
where diagonal and non-diagonal elements of the matrix are defined as 

G,,, (;, A;, r) = 4 Tr {Ii,Jk,z) Pi {I, k, r) , (35) 



G+-- {l,k,T) = {G-'+ {l,k,T))* = 4Tr (/+/-) PI {l,k,T) . 
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Figure 2: A The r-dependcnce of B.± (Afefc,r) (r = tuJE, uje/2tt - 10"Hz), for 6^ = 1/4, s = 10"^ and 
A5fc = 6^ — 6fe = — 3 • lO^'^, 1 • lO^-^, 3 • 10~^. The case B corresponds to the enlarged values of R± (Afe/^r) 
only for Afo^ = 1 • 10"^ 3 • 10"^ in the range of 10"^ < r < 10"^. 
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In the interaction representation relative Hamiltonian 



ho {I, k) = hs- {{coi (l) - a/2) h (0 + (W7 (fc) - a/2) h {k)) 
the reduced density matrix has the form 

Pin {I, k,) = exp (iha (l, k) t) p (l, k, t) exp (-i/io k) r) . 



(36) 



(37) 



Let there be the pure triplet entangled state of two removed spins I and k with zero total z-projection / = 1, 
M = 0, which belongs to the same sublattice, realized by the certain external action in the initial moment 
r = 0. It will be described by reducer state vector |1, 0) = ^JlJ2{\ t4) + | 4-t)) or following density matrix: 



PI (Z,fc,0) = |l,0)(l,0| = - 




Olio 
110 




1/4{1 - 4 {Ii^Jk,.) + 2 [irlt) + 2 [Itlk)}^ (38) 



where Gz,z{l,k,0) = —1, (j+'~(Z, fc, 0) = G~'+(Z,fc,0) = 2. The concurrence of this entangled state has 
maximum value C = 1. 

The evolution of such two-qubit state is due to the qubit interaction with magnons which in considered 
model is described by the following perturbation Hamiltonian in interaction representation 



Ahis [l, k, t) = Ahis [l, t) + Ahis [k, r) 



(39) 



= exp {iho {I, k) r) {Ahis (0 + ^^is [k)) exp {-iho {I, k)) r, 

where values Ahis {^jt)^ ^his {k,T) are determined above by expression (8). 

Let us suppose now (as previously in Scc.3), that the interaction of nuclear spins in ground coherent state 
(38) with electron system of the antiferromagnet is turning on at initial moment r = when non-disturbed 
density matrix has the form of direct product: p (0) = pi {I, k, 0) ps (0) = pi {I, k, 0) |0)(0|. 

We will write next the equation for density matrix in the interaction representation pin {I, k, r): 

idpin {I, k, t) /dr = [Ahis {I, k, r) , pin {I, k, r)] . (40) 
Prom Eq.(40) for density matrix in the second theory of perturbation theory it is follow the equation 

idpin{l,k,T) /dT =^ [Ahis il,k,T) , p,n{l,k,0)\0){0\]- (41) 

T 

-i J [Ahis{l,k,T), [Ahis {l,k,T') , PI {l,k,0)[0){0[]]dT' . 


By using Eq.(41) and making the cyclic permutation, we will obtain for the elements of reduced density 
matrix (34) the equations: 



d (G,,, {I, k, r)) /dT = 4Tri {Ii,Jk,z) dpin {1, k, r) /Or : 



(42) 



and 



-4Tr/ J {0\[[{Ii,Jk,z),Ahis{l,k,T)], Ahis{l,k,T')]piil,k,0) \0)dT', 



dG+'- {I, k, t) /dT = 4 Tr (J+4-) d Pin {I, k, r) /Or « 

r 

-4 IV y (0| [ [(/,+4-) , Ahis {I, k, r)] , Ahis {I, k, r')] Pi {I, k, 0) |0) dr'. 



(43) 



Let us take into account Eq.(38), define the tensor longitudinal component in the form 

G,,, {I, k,T) = - exp [-711 {I, k, r)] (44) 
and the tensor transverse component in the form 

G+'- {l,k,T)^2 exp [^ {cui (l) -oji{k))T- {I, k, r)] . (45) 

To obtain the expression for the rates of pair relaxation of longitudinal and transverse components in the 
context of second order of perturbation theory, we will write: 

aC,,, {I, k, t) /Ot = d7|| {I, k, r) /dT = d^l [l, k, r) /dr = (46) 

T 

= -4 TV, J (0| [[{Ii,Jk,z),Ahis {l,k,T)],Ahis {l,k,T')] (- {Ii,Jk,,) + 1/2 (l+Ij:) + 1/2 |0)dr' = 



= d7|| (Z,t) /dr + (i7|| (fc, t) /dr + dry (fc, I — k,T) /dr 

and 

Re dG+'- {l,k,T) /dr = -2Re djjL {I, k, r) /dr = (47) 

T 

= -4ReTr, J (0| Ij:) , Ahjs {l,k,T)] ,Ahjs {l,k,T')] (- {IiM + 1/2 (/+/,) + 1/2 (1^1+)) |0)dr' = 



= -2 Re d'y± {I, r) /dr - 2Red-f± {k,T) /dr - 2Redr± {kj - k,T) /dr, 
where the correlation part of longitudinal two spin relaxation rate 

df|| {k,l- k,T) /dr = 

r 

= -^"^1 j (0| [[h.,^his {I, r)] h,., Ahis {k, r')] ((/+/,-) + {lrl+)) |0) dr' = (48) 



^4Rey (0| [5+(r'),5rW] lO)'^^'- 





Consider next the correlation part of transverse two spin decoherence rate 

Rerff_L {k,l-k,T) /dT = 



T 

= -2ReTr/ J (0| , Ahis {1,t)] 4", Ahis{k,T')] {Ii,Jk,z) + 



+ [[/fe, A/izs (fc,r)] /+, Ahis {1,t')] {Ii,Jk,.) |0)dT' = (49) 

T 

= -^2 Re y (0| [S- (r') , S+ (r)] + [5+ (r') , 5," (r)] |0)dr' = 



T 

= -^4ReJ{0\ [S+{r'),Si{r)] |0)dr', 
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Figure 3: The r-dependence of decoherence rate R± {Abk, I — k.r) for values Afo^ — 3.10 '^j s = 10 ^, I — k = 
200,299,300,301. 

that is the correlation part of decoherence rate is equal to the correlation part of longitudinal relaxation rate. 

Let us rewrite next the expression (49) similarly as it would made in Sec. 3, accounting only low magnon 
excitation mode with energy : 

Redf_L {k,l-k,T) /dr^ 



(Any 



r2Re 



u* {q±, E) u {q^, Qy, E) exp (— (i£'_ — s) (r' — t)) exp [i {q^k — q'^l)] dEdq^dq^dr' 



Vl + + ^ il^J I - {E{q^) - b(^i+k)/2) sin (q^ (l - k)) + s cos (g^ (/ - fc)) 



(50) 



[{E (q±) - &(;+fc)/2) sin {{E (g±) ~ fe(i+fc)/2) t - q^{l - k)) - s cos {{E {q±_) - ^(;+fc)/2) t - q^ {I - fc))] exp (- 

{E (q^)-b^i+k)/2y 

Taking next in to account that J sin [q±_ cos 'p{l — k)]d(p — and introducing again the variable £, ^ E (g^) - 





be, we will write E {qj 



Hl+k)/2 

27r 



E {qj_) - bk - g {I - k) /2 ^ i + Abk - g {1 ~ k) /2 and obtain (Fig. 3) 



— Re dFj. (fc, ; - fc, r) /dr « i?^ (A6fc, Z - fc, r) 
3a^ 



(51) 



l + bl + bc + nY{^ + Abk-gil'k)/2,T)Jo 
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The decoherence rates of entangled qubit pair are due to decoherence of one spin states I, k and also to 
the correlation between nuclear spins I, fc. The initial diagonal and non-diagonal elements of density matrix 
1 + {I, fc, r) and G"*"'" (/, fc, r) are decreased with full rates RedFn (/, fc, r) /dr and Re(irj_ (/, fc, r) /dr. Note 
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that the asymptotic vahic for correlation part of dccohcrcnce rate (Fig. 3) may periodicaUy change sign, if 
— g (I — k) /2 < 0, that is for states after "turning point". Likewise, the indirect interspin interaction is 
described by precisely the same oscillating sign-changing function after "turning point" (see Eq. (74) in Ref.[4]). 
The concurrence for entangled two-qubit state can be obtained by using the Wootter formula (Ref.[6]) 



C {I, k, r) = 1/2 max{|G+'- {I, fc, r) | - (1 + G,,, (/, k, r)) ; 0}. 



(52) 



Taking in to account that Re7_L {k, r) = 7|| {k, r), |G+' {I, fc, r) | ~ 2exp (— Re7_L {I, k, r)), Gz^z {h k, r) w 
- exp (-7|| {I, k, r)) and also Eqs (44), (45), (47), we will obtain 

C{l,k,T) = l/2{3exp(-Re7_L {1,t) - Re-/± (fc,T) -Ref_L (fc,/- /c,t)) - 1}. (53) 
For the concurrence-damping rate we will then write 

dC {I, k, t) /dr = -3/2d (Re7_L {1,t)+ Re7_L {k, r) + Ref_L {k, I - k, r)) dr = (54) 

= -^{R± (A6z, r) + i?x {Abk, t) + {Abk, I - k, r)) . 

Note, that for A6;; > parameter Abi = Abk — 9{l ^ k) can change the sign. The expression R±_ {Abi , t) 
takes the form 

R^{Abi,T) = R^{Abk-g{l-k),T)= J (^^l + bl + bc + ^^Y{^ + Abk-g{l-k),T)d^ (55) 



Finally, for the value of concurrence damping rate we will obtain 

9a2 



where (Fig. 4) 



dC{l,k,T) /dr = -^Rf {Abk, l-k,T). 



Rl {Abk, l-k,T) = Ri_ {Abk, t) + R± (A6;, r) + Ri_ {Abk, l-k,T) 



(56) 



I 



l + b^ + bc + n- [Y{^ + Abk,T) + Y{^ + Abk-g{l-k),T) 



(57) 



+ Y{^ + Abk-g{l-k)/2,T) Jo 12 [{be + if -bl {I - k) 
The concurrence damping rate tends at r ^ oo to positive value: 



d^. 



Rl{Abk,l-k,oo) =s J {\/l + bl + bc+i 



{^ + Abkf {^ + Abk-g{l-k)f 



+ 



+- 



Jn Jl2 



{bc + O^-bl {l-k) 



{^ + Abk-g{l-k)/2Y 



di > 0. 



(58) 
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Figure 4: The r-dependence of {AbkJ - k,T) for values Abk = 3.10 ^, s = 10 ^, I - k = 
200, 299, 300, 301. 

5 Adiabatic decoherence caused by interaction with nuclear spins 
of random distributed isotopes substituting the basic isotopes in 
antiferromagnetic structure. 

Between other mechanisms of decoherence it should be pointed out the adiabatic mechanism that is determined 
by magnetic interaction of nuclear spins-qubits with electron and nuclear spins of impurity atoms, which play 
here role of an environment (Ref.[5], §5.4). 

Interaction of nuclear spins with magnetic moments of impurity paramagnetic atoms is of no concern as 
compared to the interaction of nuclear spins with electron spins of own atoms. This mechanism is largely 
suppressed for a high degree of electron spin polarization (at B/T > 30T/K). 

Let us consider here a decoherence model, where nuclear spin-qubit interact with nuclear magnetic moments 
of randomly distributed impurity isotopes in basic nuclear spin- free antiferromagnet. 

Hamiltonian of dipole-dipole magnetic interaction of considered nuclear spins has the following form: 

JV 

Hl,h^p = ^ X! ^'^P Ialp,imp {rt, t) , (59) 

where 

Do.fi in) — ^3 ^5 ^ j , (60) 

7//27r, ^i.imp/'^''^ are giromagnetic ratio of quantum register nuclear spin-qubit and of impurity isotope nuclear 
spin, Ti is radius- vector of distance from the position of nuclear spin-qubit to position of i-th impurity nuclear 
spin, ^ I ^imp^I f) ,imp (ji^t) IS fluctuating magnetic moment of impurity isotope, which produce the random local 
field: 

N 

ABa (t) = - ^Da^ (r,) {Ip^imp {ri,t) - {Ip^imp {ri))) hi- (61) 

i,/3 
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iV 

The mean value {"fiBa {t)) = imp {u)) determines the shift of qubit resonance frequency. 

The correlation function for random modulation of nuclear spin resonance frequency is determined by the 
following expression 

(Aw (r) Aw (0)) = Ci^imp / {'^h,^mp (r, r) Ifi^.mp {r, 0)) - {Ifj^imp {r)f)dr, (62) 

where Cj^imp is concentration of impurity isotopes . 

The expression for correlation function in considered case of adiabatic decoherence takes the form 

(Aw (r) Aw (0)) = (Aw2) exp {-t/T^Ump) , ■ (63) 

where quadratic mean value of modulation frequency 

(Aw^) = C/,i„p (^^ini,imph) ^0 {{Iz,imp) - {Iz,imp)^) , (64) 

{Ihmp) - {Iz,imp? = (1 - tanh^ {ji^impBn/2kTi)) /4, (65) 

a is minimal distance to impurity nuclear spin, with is of order of lattice constant ~ 1 nm. 

If it is believed that T'||,imp ~ lO'^ sec ^ Tjj (the so-called condition of rigid lattice) and that Tq > 1 sec, 
for the determination of allowable isotope concentration we will obtain the condition 

l/Tl « (Aw2) = Ci^i^p (^!^jjjj^,^^n) ' ^ (1 - tanh^ {\ji^irnp\Bn/2kTi)) < 1 sec'^. (66) 

For values B/T > 30 T/K we will obtain for allowable concentration of impiuity isotopes the highly rigid 
condition Cismp < 10^'"' cm^'', (^ 10^'"' %). However, if temperature T/ for nuclear spins of impurity isotopes 
corresponds to value, for which almost full nuclear spin polarization takes place \^i,imp\B?i/kTi > 1, that is 
T/ < 1 mK, so the allowable concentration in isotope-pure antiferromagnet will Ci^imp% < 4.5 • 10~^%. It will 
rapidly increase on further lowering of nuclear spin temperature. 

Eventually the suppressing of this mechanism calls for appropriate cleaning of substrate from impurity atoms 
and using very low spin temperatures for nuclear spins. 

6 The encoded DFS (Decoherence-Free Subspaces) logical qubits 
are constructed on clusters of the four-physical qubits 

Let us consider here two encoded DFS logical qubits \0l) and |1l) with zero total angular momentum J = 0, 
mj = 0, which states are constructed on cluster of four states of physical qubits: 

|0i) = l/2[(|01)-|10))®(|01)-|10))] (67) 

|U) = I/VS [|11) ^ |00) + |00) ^ 111) - 1/2 (|01) + |10)) ® (|01) + |10))] . 
They form two dimensional subspace for quantum operation. 

As is shown in Refs. [7-9] this subspace represent, as it is called, the strong collective decoherence free 
subspace (DFS), if the states of four physical qubits are the eigenstates of Hamiltonian with antiferromagnetic 
interaction (A > 0) of XXX type in the absent of external field 

ht,, = A/2 ^ / (fc) / (0 = A/2{^ / (A) ^ / (0 - 31}, (68) 

k^l k I 

where 1 is unit four-dimensional matrix (here symbol (8 is omitted) and is not relevant. It has eigenvalues 
A/2 J (J + 1) with J = 0, 1, 2 The DFS states (67) are the lowest state which identical to ground state of 
Hamiltonian (68) with J = 0, m = and has two-fold degeneracy. 
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We will use now the before obtained expression for effective Hamiltonian of two nuclear spins, belonging to 
common sublattice in quantum register ([4], Eq.(A2.8)) with interaction of XXO or two-dimensional isotropic 
type: 

/i..o = - E {[^^ (j') - a/2 - m h U) + U {j, j) /2- (69) 

-u {k, I) [I- (fc) /+ (0 + /+ (fc) /- (0] } 

We write next the total nuclear spin Hamiltonian of quantum register neglecting the terms W (j) « U = 
Const, of order of and take in to account that the difference of w/ (fc) — w/ (l) = ^i/^s9 {k — I) <^ a ^ 10~^ 
for (fc-0 < 1/5 ^ 10^): 



hii = - (w/ - a/2) ^ 4 (fc) - ^ [/ (fc, I) [I- (fc) /+ {I) + 1+ (fc) I- (l)] . 

k k^l 



(70) 



Our Hamiltonian (70), written for four qubits, differs from Eq.(68) in that it has identical values of interaction 

4 

parameter A for different qubit pairs and it is free from the term ^ 7^ (fc) 1^ {I). In addition, the Hamiltonian 

ki^l 

(70) has operator ^ (fc). 

fc 

In Ref.[4] it was obtained that the dependence of indirect interspin interaction U {k,l) from distance for 
(Z — fc) > 2Abk/g = 2 {be — 6 — gk) /g takes oscillating character with quasi-period which gradually decreased 
with increasing of the distant {I — k). In this case indirect interaction periodically change the sign. Consequently, 
it is possible to choose the nuclear spin position in such a way that the value of interaction U (fc, I) = —U < 
would be the same for all four qubit position in the nuclear spin chain considered as a quantum register. 



{LJi-a/2)J2lz{k)+U{l/2 



fc=0 



^ I- (fc) 1^ (0 + E (^) E i~ (0 



Lfc=0 
4 



fc=0 



21} 



(c., - a/2) h (fc) + C/{E ^ (fc) E ^ (0 - E (fc) E (0 - 21} 



(71) 



fc=0 



fc=0 



fc=0 



1=0 



The identity component of h'^^^^ in Eq.(71) is not relevant here and will be next omitted. 
The eigenvalues of such Hamiltonian are 

E{J,m,j) = -\uji - a/2\mj + U {J {J +1) -rn^j) , mj = 0, ±1, .... ± J, J = 0, 1, 2. 

They are tabulated below in Table: 



(72) 





J = 


J = 1 


J = 2 


mj 


= 





2U 


6U 


mj 


= ±1 




±\i^i - a/2\+U 


±\uji - a/2\ + 5U 


mj 


= ±2 






2±\iJi - a/2\+2U 



Because that \u)i — a/2\ U ^ a?, the states of four considered physical qubits here represent approximate 

so-named weak collective decoherence free subspacc (Rcf.[9]). The ground state of Hamiltonian (71) corresponds 
to nondegenerate state with J = 2, rrij = —2 and to eigenvalue —2\(jji — a/2| -|- 2U . It can not use for two logical 
qubit encoding. 

However the all states with raj = are identical to the strong collective DPS state with J = 0, mj = which 
overlie the ground state by 2 |w/ — a/2| — 2U. It is represented by two-fold degenerated state of Eq.(67) type. 
This state at low temperatures \ijJi — a/2\ » T/ (/tw^/fce) (T ~ 1 mK) is a metastable state and, consequently, 
can be used as two DPS encoded logical qubits. 
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Conclusion 



It was considered the results of theoretical investigations of one qubit and two qubit nonadiabatic decoherence 

and longitudinal relaxation caused by interaction of nuclear spins-qubits with virtual magnon excitations in 
antiferromagnet. It turns out that the character of decoherence processes essentially depends on antiferromagnet 
anisotropy (parameter be) and on inhomogeneity of external field (parameter g). As this takes place, the 
temperature whereby the thermal magnon excitations are excluded and the two-magnon spin-lattice relaxation 
is especially suppressed, should be defined by values T <C Tc (1 — b/bc), To = h'jsBc/kB- As an example we 
have also considered decoherenc of pair qubits maximally entanglement state and have calculated the concurrence 
damping rate. 

The other mechanism of quantum state decoherence is adiabatic process of the resonance nuclear spin 
frequency modulation caused by dipole-dipole interaction with nuclear spins of impurity isotopes. The necessary 
degree of the adiabatic decoherence suppression can be obtained at spin temperature less than ImK and for the 
concentration of impurity nuclear spin containing isotopes less than 10~^%. 

Finally it was discussed for considered model of quantum register the possibility of construction of DFS- 
states. 
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